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The Einstein equations for an isotropic and homogeneous Friedmann-Robertson-Walker Universe 
in the presence of a quintessence scalar field are shown to be described in a compact way, formally 
identical to the dynamics of a relativistic particle moving on a two-dimensional spacetime. The 
correct Lagrangian for the system is presented and used to construct a spinor quantum cosmology 
theory using Breit’s prescription. The theory is supersymmetric when written in the Majorana 
representation. The spinor field components interact through a potential that correlates spacetime 
metric and the quintessence. An exact supersymmetric solution for k = 0 case is exhibited. This 
quantum cosmology model may be interpreted as a theory of interacting universes. 
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The quantization of cosmological Universes has been a 
very active field for decades. Quantum cosmology field 
aims to construct a quantum theory for the entire Uni¬ 
verse. Several possibilities have been developed (for a 
comprehensive review, we refer the reader to Refs. [II- 
i). One way to proceed is to perform a canonical quan¬ 
tization of the classical dynamical equations for the evo¬ 
lution of the Friedmann-Robertson-Walker (FRW) Uni¬ 
verse. This is achieved by replacing the momentum by a 
derivative operator with respect to the scale factor vari¬ 
able. The final equation is widely and generically known 
as the Wheeler-DeWitt equation H, Q. This equation 
has been used to study several features of a quantum 
Universe [Ml. 

In this work, we study a FRW Universe with dark 
energy (described by the quintessence field). We show 
that there exists a deep correlation between spacetime 
and the quintessence field which allows to describe the 
Universe as a relativistic point-like particle. We study 
the quantum version of this theory, but instead of using 
the Wheeler-DeWitt equation, we perform a more fun¬ 
damental kind of quantization d la Dirac. We find that 
this quantization allows us to obtain a supersymmetric 
theory for quantum cosmology. 

The evolution of the spacetime metric inter¬ 

acting with the quintessence (massless scalar) field (jy{x^), 
characterized by a potential V(^), can be studied using 
the classical Lagrangian density 


^ 


^R- A - (f),^ +V{4>) 


( 1 ) 


where g stands for the determinant of the metric g^,^, 
R is the scalar curvature and A is the cosmological con¬ 
stant. We choose = 1 (where G is the gravita¬ 


tional constant and c is the speed of light). Notice that 
the quantities related to the quintessence field appear in 
the Lagrangian density with a sign which is the opposite 
of the one used to describe scalar fields. It is worth to 
mention that the quintessence field is described by the 
same Lagrangian density as the tlaplon field [l^ which 
describes propagating torsion. The action S = f Cd^x is 
varied with respect to the metric tensor and with respect 
to (j) to get the evolution equations. Let us consider the 
simplest case of a Universe described by the line element 
for an isotropic and homogeneous FRW spacetime [l^ 


ds^ = dt^ 


a{tf 


dr^ 
1 — 


-br2(d6|2 -bsin^ 6d(f^) 


, ( 2 ) 


with the usual definitions for the scale factor a{t), and 
the curvature constant k (where k may be ±1 or 0). The 
dynamics of the FRW in the presence of quintessence is 
obtained by varying the action S, getting the two second 
order equations 

2 -+(-)'+ 4 + = ( 3 ) 

a \a J 2 


9 + H-—— — 0 , 

a d(p 


and one first order constraint 



+ V{g>) 


= 0 , 


( 4 ) 

( 5 ) 


where we introduced V{(j)) = V{(j)) — A. Eqs. (I3|) and 
(ini) correspond to Einstein equations, while Eq. (|3|) is the 
Klein-Gordon equation for the quintessence field. An¬ 
other useful equation can be obtained by manipulating 
Eqs. m and to get 


d 
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This set of equations has already been studied taking 
advantages of its symmetries [ITI - II^ . 

The Lagrangian which is usually adopted for this sys¬ 
tem is 


and the constraint 


dp do 


( 13 ) 


L = 3ad^ — Ska — 




The Fermat-like Lagrangian Lp which gives rise to all 
(7) three equations (|3]), (HI) and ([5]) is 


and it gives rise to the dynamical equations m and 
O- However, it is important to emphasize that the 
Lagrangian © does not produce the constraint equa¬ 
tion ©• This constraint is equivalent to imposing that 
the Hamiltonian H associated to L given by 0 vanishes 
identically, i.e., 


dL. dL : 

H = T~rd H-r 

da dip 


• -L = 0. 


( 8 ) 


Now we can introduce the remarkable change of variables 


P = 


3/2 

o ^ ’ 



(9) 


which recasts the Lagrangian L in a “kinetic energy mi¬ 
nus potential energy (T— V)" form (disregarding the fact 
that both T and V have the “wrong signs” in the 0 as¬ 
sociated terms) 


L^L = ^{p^-pH^)-v{p,e), ( 10 ) 


Lf 


2Vip,0) 




(14) 


where A is, in principle, an arbitrary parameter. To re¬ 
produce the relativistic equations of motion, A is defined 
by Liineburg’s parameter choice [13, [HI 



It is straightforward to prove that varying the Lagrangian 
da with respect to p and 9 one gets Eqs. ([3]), (|3]) (when 
rewritten in terms of p and 9) while the constraint ([5]) 
arises much in the same way as = 0 ap¬ 

pears in the Lagrangian description of the dynamics of 
a relativistic particle. This statement may equivalently 
be related to the (scalar) eikonal equation in geometrical 
optics which, in turn, is equivalent to the (zero energy) 
Hamilton-Jacobi equation in mechanics. 

We can go further in the description of the FRWQ 
dynamics if we rewrite the Lagrangian (1141) as 


with the general potential 

1/3 


Lp = ^V{^,9) e2« (02-^2^, 
^(P)^) = 3 k p ^ p V{9). (11) where we have introduced the new variable 


(16) 


Note that the Lagrangian L defined in (flOll . which de¬ 
scribes the evolution of a Friedmann-Robertson-Walker- 
Quintessence (FRWQ) Universe constituted by geometry 
(represented either by p or by a) and quintessence (rep¬ 
resented by 9 or Ip), shows that the Universe evolves as a 
relativistic particle moving on a two dimensional surface 
under the influence of the potential U, thus geometrically 
unifying gravity and quintessence (dark energy). In or¬ 
der to construct a Lagrangian which in fact gives rise to 
all three equations (|3]), dH) and ([S]) it is enough to recall 
that [131 Jacobi-Maupertuis and Fermat principles give 
rise to identical equations of motion in classical mechan¬ 
ics and geometrical (ray) optics except for the fact that 
Fermat principle also produces a constraint equation. It 
is worth mentioning that exactly the same results are 
reached for the case of a relativistic particle moving on 
a two dimensional conformally flat spacetime. There¬ 
fore, the description of the FRWQ system in terms of 
a Fermat-type Lagrangian is established by defining the 
relation between the potential V{p,9), the (“refractive 
index”) conformal factor 


( 12 ) 


C = Inp, 


(17) 


and 9 = d9/dX, ^ = d^/dX. Remarkably, the previous 
Lagrangian is identical to the one of a relativistic parti¬ 
cle in a two-dimensional conformally flat spacetime, with 
metric g^i, = (where is the flat spacetime met¬ 

ric) and the conformal factor is 


H = \/Ue« 



-V{9)e 




1/2 


(18) 


Thus, the Lagrangian (1131) is written as 


Lp — 


dxf^ dx’^ 


9ti.v 


dX dX 


(19) 


with x° = 9 and x^ = 

The correspondence of the spacetime FRW geome¬ 
try and the quintessence (dark energy) field at the La¬ 
grangian level (fT6|) of a relativistic particle, allows us to 
study the quantum version of this theory. Quantizing the 
above FRWQ Lagrangian will give rise to a quantum cos¬ 
mological model. To avoid problems in the quantization 


n\p,9) = -2V{p,9), 
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procedure, we study the cases where > 0, considering 
static manifolds only where there exists a family of 
spacelike surfaces always orthogonal to a timelike Killing 
vector. This implies that dsg^v = 0, or 



( 20 ) 


Thereby, the original potential V{6) is a constant, which 
is essentially equal to a cosmological constant. Under this 
assumption, the associated Noether conservation law en¬ 
sures that 0 does not change sign (see, for instance dU for 
a (f) independent potential) and, therefore, 6 may be used 
as the evolution (time) variable which is exactly what 
the variational principle and quantization procedure sug¬ 
gest. Therefore, notice that the quintessence field acts as 
a Super-time in this new description where the particle 
is moving in an “effective” two-dimensional conformally 
flat spacetime. 

Classically, it can be rigorously shown [2^ that the 
Hamiltonian for the system described by Lagrangian dUD 
is 


H = v^\/l - + ^ , (21) 

where tt is the canonical momentum, and We 

use this Hamiltonian to construct the quantum theory 
for the FRWQ system. In order to avoid factor order¬ 
ing issues, the quantum Hamiltonian operator 'H may be 
constructed from its classical analogue m as 


n = , (22) 

where p is the momentum operator 

because to -S' = —id^. The quantum equation that de¬ 
scribes the quantization of the FRWQ system is 

(24) 


where 'k is the wavefunction for the FRWQ system. In 
principle, one may ask whether there are ways to con¬ 
struct other Hamiltonian operators that differ from (|2^ , 
giving rise to quantum theories which are not equivalent 
to the one described by ihdg^ = 'H'i’ (see, for instance 


2^). This point is subtle, and the answer is affirma¬ 


tive, however the operator (1^ has the advantage that it 
reproduces results from quantum field theory in curved 
spacetimes as we will show below. 

The quantum equation iHde'i’ = 'HT is only possible 
due to the close association between the FRW geome¬ 
try and the quintessence scalar field at a Fermat-like La¬ 
grangian level. Each quantum cosmology model requires 
a closed form for the Hamiltonian (1^^ . which basically 
consists in finding a quantization procedure for which 


the square-root of the Hamiltonian operator (|22|) can be 
explicitly written. 

Solving the square-root for a Klein-Gordon spinless 
relativistic particle [2^, one can obtain the Wheeler- 
DeWitt Super-Hamiltonian formalism. This is the most 
famous description of quantum cosmology. However, we 
can choose a different way to construct the square-root 
using matrices, in a fashion similar to the quantization 
of a relativistic particle carried out by Dirac. This idea 
has been used before to get a Dirac square-root formu¬ 
lation (departing from the Wheeler-DeWitt equation) in 
a general form in the context of quantum cosmology and 
supergravity [26l - [33| . Differently, due to the close rela¬ 
tion of theory (I^Tl) to a point-like particle descrition, in 
this work we use a more direct way to proceed with the 


quantization of this system. Breit |35| showed that there 


exists a correspondence between the Dirac and the rela¬ 
tivistic pointlike particle Hamiltonians via a prescription 
of replacement of the particle velocity and the Dirac ma¬ 
trices (as well as the prescription in Schrodinger or Klein- 
Gordon theories where the energy and momentum may 
be replaced by the time and space derivatives). Thus, 
Breit’s prescription implies a classical and geometrical 
interpretation of the spin. 

Following Breit’s ideas, we can obtain the relation 
= goo + from the classical Hamiltonian (l?T|) . From 
that relation we can obtain 


^ = (25) 

where we defined the velocity variable q = -k/H. Using 
this variable, the square-root of the Hamiltonian (l2T|l can 
be written 

H = ^o{^ + V^T^) . (26) 

Breit’s interpretation [3^ corresponds in the identifica¬ 
tion of the Dirac matrices as g —>■ a, and i/l — —>■ /3. 

Breit showed the these identifications are consistent with 
the postulates of Dirac electron’s theory. The implica¬ 
tions of this prescription have been investigated with the 
purpose of understanding the underlying nature of the 
spin or antiparticles [s^. 

Following Breit’s prescription, the Hamiltonian opera¬ 
tor can be written using Dirac matrices (a and P) 
as 


= ia-p + (27) 


where a and /? are the flat spacetime Dirac matrices, as 
the curvature is already taken into account in D. Because 
of our FRQW system is two dimensional, the matrices 
are two-dimensional as well (as Dirac matrices are two- 
dimensional in 2 or 3 dimensional spacetimes). Using the 
operator (1271), the quantum mechanical equation = 
reads (with H = 1) 


I'j 




■ I'j 



119 In D \ 
2~^ ) 




(28) 
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where now 4/ is a spinor, 7 ° = /3 and 7 ^ = 7 °a. Because 
of the dimensionality of the FRWQ system, there exist 
only two 2 x 2 matrices satisfying the algebra { 7 ^, 7 '^} = 
277^“^. 

The strongest feature of Breit’s ideas and quantization 
procedure (|T7l) is that the Dirac equation (|2^ is exactly 
the same equation obtained from Quantum Field The¬ 
ory in curved spacetime [s^ for the Dirac equation in a 
two-dimensional conformally flat space. In general, the 
curved-space Dirac equation is 

= (29) 

where the vierbein e^a are defined through = 

(and = ^b)- The spin connection is 

defined as 




— e^ue'^b,n + , (30) 

where are the Christoffel symbols. In a two- 

dimensional conformally flat spacetime, the zweirbein 
components are eo° = D, ei^ = O, e°o = I/H, and 
e\ = I/D. Thus, Wabo7“7^ = 2woio7'’7\ Wa6i7“7^ = 
2 a;oii 7 ° 7 ^, wqio = T^io = 9jlnD, and wqh = F^n = 0. 
Using these results, we recover Eq. (1^ . The equivalence 
of these two different methods reinforces the Breit’s pro¬ 
cedure. 

It is necessary to remark that quantum cosmological 
theories in a two-dimensional Lorentzian geometry can be 
derived in other ways different than using a quintessence 
model. For example, this has been studied in a spinor 
model for two-dimensional dilaton gravity [s^, or in the 
spmor representation of the Wheeler-DeWitt equation 

Furthermore, with Eq. (l28l) . the expected value of the 
scale factor can be readily calculated as 

{am = , (31) 

evaluating the integrals and 'h in terms of a. Here, 'I'l is 
the transpose conjugated of the wavefunction 4', and the 
probability density of the Dirac field is 

pa 1 

|^i(6»)P = / (32) 

Jo a 

We can see that (a) does depend explicity on the coupling 
between wavefunctions, and also on the Super-time 9. 

Now, let us define the new wavefunction 4> = 

•\/D4> exp)^^^), which satishes 

(-f7^ ^ -b d) $ = E;7°$ , (33) 

which can be found using Eq. (1^ . Eq. (1551) couples the 
components of the spinor $. Defining the operators P+ 
and P- , and the spinor $ as 



P+ = 




then P|4> = (p+, and = ip-. The system (1551) can 
now be recasted in the two coupled equations 


Q±p± + W±p^=0, (35) 

with the operators 

Q± = -^{PWP±)-^ + n-E{plJ°p±), 

W± = -^(plJ^P^)^ - E(pIj°P:^) . (36) 

Next, let us focus in the different choices for the matri¬ 
ces which satisfy { 7 ^, 7 “"} = 2^^", giving different types 
of evolution equations. One interesting posibility is the 
choice of 

Pi7'PT=0. (37) 


In this case, it will be fixed the form of the matrices to 
be 


7 


0 




(38) 


also impliying that 

Pl7°P± = 0 , Pl7^P± = Tf, Ph^’PT = 1 ■ 

(39) 

Matrices (I38|) correspond to the two-dimensional Ma- 
jorana representation, which implies that $ is real. Be¬ 
sides, the operators (1551) reduce to 


Q± = T:^ + ^, W± = -E. 
Thus, dynamical equations (1551) become 


(40) 


(t^ + n'J<p±= Ep^ . (41) 

Eqs. (15T1) correspond a set of supersymmetric equations 
of quantum mechanics ii l4ll - l^ . This Supersymmet¬ 
ric Quantum Cosmology (SSQC) theory can only be ob¬ 
tained in the Majorana picture. This SSQC scheme is 
different from previous models d, d, , as here the 

SSQC emerges naturally owing to the correspondence of 
the FRW spacetime and the quintessence. As Eq. dSD 
has a supersymmetric structure, the two spinor compo¬ 
nents are the super-partner of each other. Each wave- 
function satisfies 


H±p± = E\± , (42) 

where the Hamiltonians operators and potentials are 

H± = -d] -b W± , W± = ±d{D -b , (43) 

respectively. As usual [i^, we can define the Super- 
Hamiltonian and the supercharge operators as 


( H+ 0 
\ 0 H- 


Q 


0 Q+ 
0 0 




0 0 \ 

Q_ 0 j ’ 

(44) 


F_ = 


(34) 


H = 
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repectively, with operators defined in Eq. dMl). Thus, we 
can find 

. 

The supercharge operators change bosonic (fermionic) 
states into fermionic (bosonic) ones. The above oper¬ 
ators have the algebra 

H = {Q, Q+} , ^ [H, Q] = 0 = (Qt)2 = p, Qt]. 

(46) 

The Hamiltonians H± has the same bound state spectra 
= E+ (for n = 0,1,...), and the bound states eigen¬ 
functions of the Super-Hamiltonian are related simply by 
and 

Besides, the energy E^ = 0 of the ground state of i/_, 
has the wavefunction 

(f) = exp • (47) 

The Majorana SSQC theory described by Eqs. dUD is 
general. An exact solution for the flat curvature k = 0 
case can be obtained explicitly. In this case 

n = . («) 

provided that we can choose H < 0. This gives rise to 
the Morse potential 

>V± = ±2e"^"^ -f , (49) 

with the change of variable x = — £ + £n and V = 
— (8/3) exp(—4^o)- It is well known (S^ . jsl ] that the 
Morse potential is one of the solvable “shape invari¬ 
ant” potentials in Supersymmetry. This kind of po¬ 
tentials fulfills yv+{x,e) — yV-{x,ei) = C{ei), where 
£r is a set of parameters, Si = /(e), and C(ei) is 
a constant. For the shape invariant potentials the 
bound state spectrum of i7_ is completely determined 
in terms of the sum of constants E~ = 

(s^ [H^. In this form, the n-state wavefunction associ¬ 
ated to iJ- is completely determined [H] by = 

and hence the 

n-state wavefunction associated to i7+ can be found. 

We find an explicit solution for fc = 0, where the so¬ 
lution (j47p for the wavefunction with vanishing energy 
becomes simply 

/_K)=exp('-iy7|E,«V (50) 

Notice that this expression represents a bound state solu¬ 
tion because <(3(!_(^ —>■ oo) = 0 (where ^ oo is equivalent 
to a —oo). Also, in general, Eq. (l42|) for </?_, reads 

+ 2e-2" - e-4") = 0 , (51) 


in terms of x. The solutions for the Morse potential are 
known. They represent a diatomic molecule with a min¬ 
imum equilibrium distance bound equal to —^o- It is 
straightforward to prove that in the region — oo < ^ < oo 
(or 0 < a < oo), the only possible solution for bound- 
states is for A = 0 [13, . 

On the other hand, there are no bound-states solutions 
for (p+. This is due to the fact that the potential W+ = 
2e“^“ -|- which appears in Eq. (1421) for (/?+ is positive 
everywhere. 

For other curvature cases /c = ±1 a deeper analysis is 
required and this task is left for future works. Neverthe¬ 
less, the SSQC theory (1^ is interesting as its origin is 
the form of the Lagrangian (TlB)) . which couples gravity 
and quintessence (dark energy) field. It is for this reason 
that the Majorana SSQC evolves in 9 and ^ (and not in t 
and a as other models in SSQC). As far as we know, the 
association between variables made in (1161) has not been 
envisaged previously, producing a SSQC theory different 
from the Wheeler-DeWitt ones. As an aside, we would 
like to mention that other Quantum Cosmology models 
are possible in this theory. For example, by choosing 

Pi7'P±=0, (52) 

the matrices are fixed to be in the Dirac representation, 
implying that 

p|7°72^ = 0, pIj°P±=TI, pI^^P^ = ±1. (53) 
In this representation, the operators become 

Q± = n±E, W±=Ti^. (54) 

Hence, it is straightforward to see that the Dirac repre¬ 
sentation does not produce a SSQC theory. 

Finally, we stress that the Dirac-like equation (l28l) 
describes a spinor quantum cosmology, representing a 
generalization of the Wheeler-DeWitt formalisms. The 
spinor components are coupled, and in comparison with 
Dirac theory for particles, we can argue that the spinor 
components represent two interacting Universes, repre¬ 
sented by the components ip±. Thus, the spinor could 
work as a description for a Multiverse. The interaction 
between the wavefunction of the Universes gives the dy¬ 
namical evolution of the expected value of the scale fac¬ 
tor. 

For example, the physical meaning of the solution for 
k = 0 shown above can be envisaged as a Multiverse the¬ 
ory of Universes behaving as a diatomic molecule in a 
bound state with vanishing energy and a minimum equi¬ 
librium distance directly related to the quintessence po¬ 
tential by (1/4) ln(—3U/8). 

This Multiverse theory is different from the Everett’s 
Many-worlds interpretation of quantum mechanics 
1^ . as each spinor component could be interpreted us¬ 
ing Everett’s theory. The Many-worlds interpretation of 
Multiverses have been largely studied (the literature is 
vast and some few examples are in Refs. [bSl - f^ b 
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The Multiverse interpretation of the wavefunction in 
Eq. (I28|) is not new in quantum cosmology (see for ex¬ 
ample Refs. However, unlike other theories, 

the mathematical equivalence between this SSQC the¬ 
ory [Eqs. (l28l) or (l33)) ] and the Dirac theory for spin rel¬ 
ativistic particles, make the Multiverse interpretation a 
compelling issue. 
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